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Abstract 

A consistent device model to describe current- voltage characteristics of metal/insulator /metal 
systems is developed. In this model the insulator and the metal electrodes are described within 
the same theoretical framework by using density of states distributions. This approach leads to 
differential equations for the electric field which have to be solved in a self consistent manner 
by considering the continuity of the electric displacement and the electrochemical potential in 
the complete system. The model is capable of describing the current-voltage characteristics of 
the metal/insulator /metal system in forward and reverse bias for arbitrary values of the metal/ 
insulator injection barriers. In the case of high injection barriers, approximations are provided 
offering a tool for comparison with experiments. Numerical calculations are performed exemplary 
using a simplified model of an organic semiconductor. 
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I. INTRODUCTION 



Injection is an important factor in many electronic devices, especially in applications 
with insulators, where virtually all charge carriers have to be injected from the electrodes. 
Examples that have been the subject of interest are electronic devices built with organic 
semiconductors such as organic-light-emitting-diodes (OLEDs). Organic semiconductors 
show many properties of dielectric materials like relatively large band-gaps or the absence 
of intrinsic charge carriers. For the processes describing the conduction in organic semicon- 
ductors itself, sophisticated models were recently developed lOj. The 
proposed injection mechanisms are taking into account the stochastic, hopping character of 
the transport, the field-dependence of the mobility of charge carriers, the roughness- induced 
energetic disorder at the interface and the specific density of states (DOS) characteristic of 
these materials. However, the space charge effects were accounted so far only concerning 
the transport inside the organic materials while the description of injection itself was based 
on the classic works of Richardson and Fowler- Nor dheim III, [12] considering injection as 
a single electron process. A comprehensive numerical model was developed by Tutis et al. 
[13I comprising the hopping transport in single- and bilayer devices and tunnel injection 
from electrodes accounting for space-charge effect. However, besides the injection barriers 
given by the bare difference of the chemical potential in the metal and the lowest occupied 
orbital (LUMO) in organic material this model involves a range of artifacts like tunneling 
factors, effective attempt frequencies etc. Being a very sofisticated numerical tool and pre- 
senting good agreement with experiments this approach does not allow analytical fitting of 
current- volt age characteristics which gives insight in major mechanisms controlling injection. 

In this work we develop a device model for a metal/insulator/metal system including 
the description of charge carrier injection at the metal/insulator interfaces. This injection 
model takes into account the electrostatic potential generated by the charges injected. The 
treatment of this problem faces the problem of self consistency, since the amount of injected 
charges depends on the height of the injection barrier, while the electrostatic potential gener- 
ated by this charge modifies the height of the injection barrier itself. This problem is solved 
by defining the boundary conditions far away from the interface, where the influence of the 
interface can be ignored. As a result, charge carrier densities and electric field distributions 
in the respective medias are coupled and depend on the conditions of the system. By de- 
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scribing the electrochemical potential as well as the dielectric displacement continuously all 
over the entire system the electric field and charge-carrier distributions can be calculated. 
One can model current- voltage- (IV) characteristics in which either charge injection or the 
transport through the insulating layer determines the current. The model includes barrier 
lowering, but this lowering differs significantly from the lowering due to the image charge 
potential which is a direct consequence of the one electron picture. 

The problem considered in this paper includes the presence of both, injecting and ejecting 
electrode with an organic semiconductor in between and, hence, a built-in potential in the 
system. Both interfaces are considered self-consistently involving no additional characteris- 
tics but the bare injection barriers given by the difference of the chemical potentials in the 
metals and the LUMO. Electric field distributions in the device are calculated analytically 
for different barriers heights in equilibrium. Then the field distributions are calculated nu- 
merically in the presence of a steady-state current. It is shown that these field distributions 
are qualitatively different for injection barriers smaller and larger than some characteristic 
value. Using the known solution for the electric field at a fixed current the current-voltage 
characteristics are computed. They exhibit, depending on the barrier heights and bias, dis- 
tinct areas of linear, exponential and quadratic dependence. In the case of devices with 
large barriers an analytic approximation for the field distribution is derived which describes 
the numerical results with high accuracy. The obtained approximation is similar to the 
well-known Mott-barrier formula but includes an effective barrier height which results from 
the consistent treatment of the metal/dielectric interface. 

II. THE MODEL 

Let us consider an insulator of thickness L sandwiched in between two metal electrodes. 
The insulator is supposed to be extended over the space with — L/2 < x < L/2, whereas 
the metal electrodes are extended over the semi spaces with x < —L/2 and x > L/2, 
respectively. In the following sections the theoretical model describing the insulator and 
the metal electrodes is introduced. The description is based on the specific density of states 
(DOS) distributions of the employed media. Differential equations will be derived, describing 
the electric field distribution by the requirement of local chemical equilibrium. 
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A. The Electrodes 

Assuming the free electron model and the Thomas- Fermi approximation [ijj], the elec- 
trochemical potentials of the metals as a function of the spatial coordinate x read, 

h 2 

«m(*) = ^l(3vr 2 n± (x)) 2 / 3 - e<f>(x) + £± (1) 

Here are the electron densities, e is the elementary charge and m ± the effective electron 

masses in the metals, (f>(x) the electrostatic potential and the bottom of the conduction 

bands. All quantities with indices ± are assigned to the electrode on the right (x > L/2) 

and the left hand side (x < —L/2), respectively. The steady state current density j is given 

by the conductivity a and the derivative of k(x) 15]. Since we consider a one-dimensional, 

monopolar case, the current remains constant across the whole space, 

cr dn(x) ,„ s 
j = i- 1 = const. (2) 

e ax 

The conductivities of the metals are proportional to the electron mobilities /i^j and the 
charge- carrier densities far away from the contacts, n^: = e/i^jn* . Due to the charge 
transfer between the metals and the insulator, space charge regions emerge, which modify 
the electric fields F^(x) also in the metal. According to Gauss's law, the derivatives of the 
electric field are proportional to the excess electron densities 5n ± (x) = n^{x) — n^: 

F*{x) = --&n ± {x). (3) 

where eo is the permittivity of free space. 

Since electron densities in metals are rather high, the value for the excess electron densities 
is small in comparison with the background electron densities, |5n <C n^. Assuming 
the linear Thomas-Fermi approximation a straightforward calculation leads to differential 
equations for F^, which read: 



j 2e /c ± 



K oo T?±"/„\ , Z7± 



x)+F±(x). (4) 



Here are the chemical potentials in the metals at an infinite distance from the 
metal/insulator interfaces. There space charge effects vanish and hence gradients of F^(x). 
With this boundary condition the solutions for F^(x) read, 



with 

being the Thomas- Fermi lengths, defining the typical length scales of the metals (Itf — 
10~ 10 m). It is clear from expressions §5§ that our assumption of the metal semi-spaces is, in 
fact, not necessary. With the very short length Itf the results are valid for any experimental 
metal thickness. In Eq. the electric fields in the metals at the metal/insulator interfaces 
F^(±L/2) remain as the only unknown quantities. 



B. The Insulator 



The energetic difference between the electrochemical potentials in the metal electrodes, 
at an infinite distance from the contact, and the bottom of the conduction band in the 
insulator is defined as the injection barrier A 1 * 1 . Accordingly the energetic difference of the 
A ± s is related to the difference between the electrode work functions E^, 

E A -/\- = E+-/\+. (7) 

The introduction of a bottom of the conduction band in the insulator means that one can 
calculate charge carrier densities by using Boltzmann statistics. 

/x 7 /-r-, a x f k s (x) + ed>(x) — E\ 
n s (x) = j g s (E - A - «" ) exp f ; ^ } j dE, (8) 

— oo 

Here, g s (E) is the density of states of the insulator with g s (E < 0) = 0, T is the absolute 
temperature and k is the Boltzmann constant. The energy scale, E, has been adjusted to 
the bottom of the conduction band in the left metal electrode. Thus, the electrochemical 
potential k s can be expressed in terms of the charge-carrier density n s , 



'n s (x) 

K, s yx) = Kl in I 

where Af is defined by 



k s (x) = kT In I j + A- + - e<f>(x), (9) 



M= J g s (E)exp(-E/kT)dE (10) 

— oo 

and can be understood as the effective total density of states available in the insulator at a 
given temperature T. One should realize that the T dependence of TV becomes weak in the 
case of a narrow-band material. 
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For wide band gap insulators all charge-carriers in the extended states, n s , are excess 
charge- carriers appearing in Eq.([9]). They have to be considered in Gauss's law, 

F^x) = -—n s (x), (11) 

where e is the relative permittivity of the insulator. With Eqs.(j2]), ([9j), f TTTT) and a s (x) = 
efi s n s {x) a differential equation for the electric field in the insulator, F s , is obtained, 

i kT 

— = F' s '(x)-F s (x)F:(x), (12) 

fi s ee e 

where fi s and a s are the electron mobility and conductivity in the insulator, respectively. 
Alternatively, this differential equation could be derived applying the drift-diffusion model 



and the 

systems 



instein relation, which relates the mobility and the diffusivity in non-degenerate 







C. Boundary conditions and self-consistency 

The self consistent treatment of charge transport through an insulator sandwiched be- 
tween two metals requires continuity of the electrical displacement and of the electrochemical 
potential across the whole system. In particular one may write: 

k(x) = continuous, (13) 
eF(x) = continuous. (14) 

These conditions have to be fulfilled particularly at the metal/insulator interfaces, eliminat- 
ing the unknown integration constants. It follows from Eq. (TT4"]) that the electric fields in the 
metals and the insulator at the interface are related by: F s {±L/2) = F±{±L/2)/e. From 
Eq. fll3l) nontrivial boundary conditions for Eq. (jl2j) are obtained: 

+ £ T #^/2) = T-^i, (15) 

which depend on parameters of both the insulator and the metal electrodes. 

Indeed, these boundary conditions are virtually independent of the net current density. 
In Eqs. flToT) the current density j is multiplied by small factors l^ F / fx^n^kT and hence, the 
boundary conditions are independent of j in most practical cases. Under this approximation 
the density of injected electrons at the interfaces reads, 



n 



: (±L/2) = Wexp 



-A±/kT±e 6 -^F s (±L/2) 



(16) 
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Hence, the dependence of F s (±L/2) or n s (±L/2) on the current is mainly due to Eq. (fl2l . 
From Eq. ffl6l) it is evident that the height of the injection barrier is modified by the electric 
field at the interface. This leads to the definition of an effective injection barrier, 

A% } = A ± =f eel± F F s (±L/2). (17) 

The barrier modification is a direct consequence of the charge transfer from the metal to 
the insulator and corresponds to the potential energy, electrons gain or lose in the electric 
field of the metal. 

III. PHYSICAL AND NUMERICAL ANALYSIS 

We now show the results of analytical and numerical calculations. The section is organized 
as follows: In the first part we show results for calculations in equilibrium where charge 
carriers diffuse from the metal electrodes into the insulator. In the second part we calculate 
electric field distributions in steady state and finally we analyze the resulting current- voltage 
(IV) characteristics. 

First, the material parameter of the metal electrodes and the insulator will be specified. 
As mentioned in the introduction we chose a simplified model of an organic semiconductor 
as example for an insulator. These materials are characterized by band-gap energies ranging 
form 1 to 3 eV making thermal excitation to conduction states virtually impossible. As in 
insulators, charge transport is therefore dominated by excess charge carriers. The DOS of 
organic semiconductors is believed to have a gaussian shape [nf, impeding in general the 
applicability of Eq. (!8]) and requiring the use of Fermi-statistics. However, in the case of 
weak disorder charge carrier trapping in tail states of the DOS is negligible and Boltzmann- 
statistics can be applied. We assume that the material specific parameters for the organic 
semiconductor adopt the typical values of: 

Af = 10 21 cm" 3 at T = 300 K, e = 3 and fi s = 10 _4 cm 2 /Fs. (18) 

An organic semiconductor with a thickness of L = 100 nm will be considered, typical for 
applications in OLEDs. As examples for the metal electrodes, calcium, barium and magne- 
sium are chosen. These metals are characterized by their electron densities in the conduction 
band n^, their chemical potentials at infinite distance from the contact Kqo, which at T = 
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equal their Fermi-energies, their electron mobilities /i m and their electron work functions Ea- 
The values for the material parameters are summarized in TabJH Here, the barriers A 1 * 1 are 
given by the energetic difference between the transport level in the organic semiconductor 
and of the respective metal electrode. 

TABLE I: Parameter values for the different metal electrodes considered in this work. 





noo [W 22 cm 3 ] 


Koo [ eV] 


Hm [cm 2 /Vs] 


E A [ eV] 


Ca 


2.6 


4.68 


66.7 


2.87 


Ba 


3.2 


3.65 


5.07 


2.7 


Mg 


8.6 


7.13 


16.9 


3.68 



In the following the solutions for the equilibrium situation will be discussed, where space 
charge zones are formed due to diffusive charge-carrier transport from the electrodes into the 
insulator. Then a constant external voltage will be considered, driving a constant current 
through the insulator, and the calculated IV-characteristics which are crucial for comparison 
with experiment will be presented. 

A. Equilibrium 

In equilibrium, the current density j vanishes and Eq. (fT2l) can be integrated: 

^F;(x) + iF s 2 (x) = -2^ 2 A 2 , (19) 

where A is a constant. The derivative of the electric field has to be negative due to Eq. ffTTT) . 
At some position between the contacts the electric field vanishes, thus, A 2 > and an 
analytic solution for the electric field can be found, 

2kT 

F s (x) = —A tan (Ax + A) , (20) 

with A being the integration constant. This solution describes the following processes: charge 
carriers are injected from both metal electrodes and diffuse far into the bulk of the organic 
semiconductor. The resulting space charge generates an electric field compensating - in 
equilibrium - the diffusion current. While the widths of the space charge zones in the metals 

8 



are of the order of Itf, they can be extended over the whole insulating layer. This means, 
that at some point the diffusion currents caused by electrons injected from both electrodes 
compensate each other: here the electric field is zero. 

From the analytical solution of Eq. (|19|) . the boundary conditions (|15[) can be simplified 
to a nonlinear system of equations determining the integration constants A and A: 



In 



'ee 2kTA 2 



+ — ± 2eAl TF tan(±AL/2 + A) = (21) 



J\f e 2 cos 2 (±AL/2 + A) J kT 

In general, Eqs. (l21l) have to be solved numerically. However, for the situation of symmetric 
electrodes the solution for the electric field is antisymmetric with respect to x — 0, which 
means A = 0. In FigJT]we show the solution for the electric field of an organic semiconductor 
sandwiched between two calcium electrodes. The injection is supposed to be barrier free 
(A~ = A + = 0). The introduction of nonvanishing injection barriers leads to a reduction 
of the electric field at the interface, although the dependence of the electric field on the 
coordinate is still given by Eq. (120]) . 
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FIG. 1: Distribution of the electric field F in units of Fq = kT/el TF as a function of the coordinate 
x in units of l^ F in case of equilibrium (j = 0) and barrier free injection (A~ = A + = 0) from, 
two Ca electrodes. Notice that the x-axis is extended between —L/2l TF and +L/2l TF ~ 610 and 
interupted between x = —490 and x = +490. 



Now the electrode at x = +L/2 will be changed from calcium to barium, to demonstrate 
the consequences of an unbalanced charge- carrier injection. The injection process at the 
barium side of the device at x = +L/2 is supposed to be barrier free (A + = 0), the injection 
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barrier for the calcium/organic interface at x = —L/2 is then calculated with the help of 
Eq.([7j) and reaches the value of A - = 0.17 eV. The effect of an unbalanced charge-carrier 
injection can be seen in Figf2J A significant asymmetry in the computed distribution of the 
electric field can be observed. 
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FIG. 2: Distribution of the electric field F in units of Fq = kT/el TF as a function of the coordinate 
x in units ofl TF in equilibrium (j = 0) with an injection barrier of A~ = 0.17 eV for a Ca electrode 
at x = —L/2 and a vanishing injection barrier for the barium electrode at x = +L/2. Notice that 
the x-axis is interupted between x = —490 and x = +490. 



Integration of the electric field over the whole system results in the built-in voltage Vbi, 

oo L/2 

-V B i= J F(x)dx = F s (-L/2)el TF + J F s (x)dx + F s (L/2) e l+ F (22) 

-oo -L/2 

using Eqs.(j5J) and (JHJ). It looks as if the integration runs virtually over the effective length 
of the device 

L eff = L + e(l TF + l TF ). (23) 

This length will be used later by the evaluation of the voltage on the device in a steady-state. 
The integral in Eq. (122)1 can be calculated analytically using the solution (120)) and boundary 
conditions (I2T1) : 

eV BI = E A -E\ (24) 

and is independent of the values for the integration constants A and A. One should be aware 
of that the built-in potential drops over the insulator and the metal electrodes. 
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B. Steady state 

Now, the steady state situation is considered, where a constant current j flows through the 
system. In this case, the differential equation (fT2"j) . describing the electric field distribution 
in the organic media has to be solved numerically with respect to the boundary conditions 
defined by Eqs. (fTBT) . 

The solution for the electric field of an organic semiconductor sandwiched between two 
calcium electrodes is shown in Figj3j The calculation was carried out for barrier free injection 




FIG. 3: Distribution of the electric field F in units of Fq = kT/el^ F in steady state for current 
densities of \j\ = 10 mA / err? (solid line) and 100 mAj cm 2 (dashed line) as a function of the coor- 
dinate x in units ofl TF . The charge injection is supposed to be barrier free with calcium electrodes 
on both interfaces. We comment that the maxima of the electric field near the interfaces are cut 
off in this plot and do not change in value with respect to the equilibrium. 



at both metal/organic interfaces (A = A + = 0) and steady state current densities of 
|j| = 10^ and 100^. In case of vanishing injection barriers many charge-carriers are 
injected into the organic semiconductor. This requires strong electric fields to compensate 
the diffusion currents near the interfaces to such an extend that the net current is j. This 
field is positive on the injecting electrode and negative on the ejecting electrode. Thus, 



18H , where the 



there has to be a position xo, often referred to as the virtual electrode 
electric field changes sign. When the current increases, this position shifts towards the 
injecting electrode. It can be seen from Figj3] that xq shifts from approximately 25 nm to 
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10 nm once \j\ is increased from 10 to 100 ^^t - In between the two electrodes, the electric 
field of the barrier free case follows approximately the electric field of a space charge limited 
system assuming ohmic boundary conditions [Uj. In this situation the dipole layer at the 
left electrode acts as a source of charge-carriers, the right electrode can be understood as a 
sink. 

The introduction of non vanishing injection barriers reduces the amount of charge- carriers 
present in the organic semiconductor and therefore the importance of diffusion. When 
the injection barrier at — L/2 exceeds some critical value no compensation of diffusion is 
necessary and the electric field near the interface becomes negative. The charge-carrier 
reservoir is exhausted and no virtual electrode is formed. In FigJH the field distribution 

0.0 
"o -0-5 
-1.0 
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FIG. 4: Distribution of the electric field F in units of Fq = kT/el^p as a function of the coordi- 
nate x in units of l^p in steady state for current densities of \j\ = 10 mA/cm 2 (solid line) and 
100 mAj cm 2 (dashed line). The organic layer of thickness L = 100 nm is contacted with calcium 
and barium with injection barriers of A~ = 0.35 eV and A + = 0.18 eV respectively. 

of a system is shown, where an organic layer with larger band gap forms a contact with 
calcium and barium electrodes. The injection barrier heights assigned to the contacts are 
A~ = 0.35 eV and A + = 0.18 eV, respectively. The electric field strength is virtually 
independent of x while its value depends on the current density. The absence of extended 
space charge zones in the insulator indicates that the system is in an injection limited mode. 
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C. Current- voltage-characteristics 



Knowledge about the distribution of the electric field allows one to determine the voltage 
drop V for a given current density j and hence to calculate the IV-characteristics. After 
solving the differential equations for the electric field for a given current density, the voltage 
applied to the system can easily be calculated by integration. For the steady-state case, 
however, one cannot define the voltage drop over the infinite range of integration as in 
Eq. (f22l) since the electric fields in the metals, Eq.Q, are asymptotically constant and do 
not vanish. To compare with experimentally observed current-voltage characteristics, one 
has to account for the built-in voltage, Eq. (1221) . which drops virtually over the effective 
device length, Eq. f[2"3"j) . and can define the voltage drop over the same length: 

L/2+el+ F 

V = - J F(x)dx - Vbi (25) 

-L/2-el- F 

In FigJS]we present IV-characteristics of a Ca/organic(100 nm)/Ca system. The injection 
barriers are varied from eV to 0.4 eV. The IV-curve calculated for vanishing injection 
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FIG. 5: IV-characteristics for an organic layer with thickness L = 100 nm contacted with calcium 
electrodes on both interfaces and different injection barriers of A~ = A + = 0; 0.2; 0.3; 0.35; 0.4 eV. 

barriers represents a limitation of the current flow through the system. At voltages up 
to ~ 0.3V the IV-characteristic shows an ohmic-like j ; ~ V dependence. We emphasize 
that this ohmic behavior is not due to intrinsic charge-carriers but a result of charge carriers 
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injected from the electrodes in the organic semiconductor even in equilibrium (Fig{TJ). Above 
~ 0.3V the system shows the Mott-Gurney j ~ V 2 dependence of the current on the voltage. 
Here the system is obviously in a space charge limited mode. The IV-characteristic does 
not change considerably up to a barrier height of A cr ^ ~ 0.27 eV. The current flow is for 
all voltages limited by the transport through the bulk of the organic semiconductor and the 
injection process has only a secondary effect. Introducing injection barriers > 0.27 eV, the 
current decreases by orders of magnitude. Nevertheless, all IV-characteristics show the same 
j ~ V ohmic-like dependence in the low voltage region like the curve calculated for barrier 
free injection. For high As this ohmic regime is followed by an exponential increase of the 
current on the voltage. For even higher voltages all curves approach the space charge limit 
and coincide with the barrier free IV-curve. 

In Figj6]we introduce a built-in voltage by changing the metal at x — +L/2 from calcium 
to magnesium. The difference in the work functions corresponds to a built-in voltage of 




vm 

FIG. 6: IV-characteristics for an organic layer with thickness L = 100 nm contacted with a calcium 
electrode at x = —L/2 and a magnesium electrode at x = +L/2. The injection barrier of the 
calcium electrode is varied from eV to 0.4 eV with an increment of 0.1 eV. The work function 
of the metals differs by an amount of 0.81 eV, accordingly the injection barrier of the magnesium 
electrode is higher by that value. 

Vbj = —0.81 eV. The injection barrier on the calcium electrode at x — —L/2 is varied from 
to 0.4 eV. Correspondingly, injection barriers on the magnesium electrode at x = +L/2 
range from 0.81 to 1.21 eV. 
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The computed IV-characteristics for high injection barriers can be understood by con- 
sidering some simple approximations. In the case of high injection barriers the amount of 
injected charges is small and the electric field in the insulator can be conveniently decom- 
posed into a constant mean value and a space charge part which is assumed to be comparably 
small, 

F a (x) = F s0 + f(x) where |/(ar)|< \F s0 \. (26) 

The mean value of the electric field is determined by the voltage drop on the effective length 
of the metal/insulator/metal system, (I25p : 

V + Vbi 



J eff 



(27) 



which implies that the integral of f(x) over the system length equals zero. 

The decomposition of the electric field represented by Eq. (j26|) leads to a linearized form 
of the differential equation for the electric field in the insulator: 



kl n s eeokl 



0. 



(28) 



which should be supplemented by the linearized boundary conditions, Eqs. ffTBT) . Taking into 
account that in the considered systems \% F jL are as small as 10~ 3 , the solution of Eq. 
reads 



f(x) = B 



2kT . / eF s0 L 

smn ( — — | — exp 



eF s0 L V 2kT 
here the quantity B is an abbreviation of 

NkT 



eF s0 
' kT '' 



J 



H s 66qF s0 



B 



2ee F s0 smh(eF s0 L/2kT) 
where we assume, approximately, 



A 



exp 



eff 



A 



kT 



exp 



eff 



kT 



(29) 



(30) 



A 



eff 



A ± =f eel^ F F s0 . 



(31) 



The analytic solution results directly in the IV-characteristic of a system with high in- 
jection barriers: 



J 



eNnsF s 



,s0 



2 sinh(eF s0 L/2A;T) 



exp 



A tf/ 



exp 



eF s0 L A 



eff 



(32) 



2kT kT J r \ 2kT kT 

This relation is symmetric in the injection barriers and is capable of describing the current 
flow in both possible directions. In the following we focus on electron transport from the 



15 



left to the right electrode. For such situations it is convenient to examine a nonsymmetric 
form of Eq. (1321) which reads: 



j - _en V + VBI Ar c:LV ( _^£l) 6Xp Lfe m) 

L *ff V kT Jexp^^i^)-!- (33) 

Equation fl33l) resembles the current- volt age characteristic for Mott-barriers [12J. An essen- 
tial difference involves the earlier introduced barrier lowering in the electrodes. In FigJT] 
we compare the approximate solution of the injection problem, Eq.( |33|) . with the numerical 
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FIG. 7: IV- characteristics for an organic layer with thickness L = 100 nm contacted with a calcium 
electrode at x = —L/2 and a magnesium electrode at x = +L/2. The stars and the triangles 
correspond to the the IV- characteristic calculated for forward and reverse bias, while the solid and 
dashed lines are the corresponding approximate solutions given by Eq. A33\). 



solution. The agreement is perfect for reverse and forward bias until the space charge effects 
become dominant and Eq. (1281) loses its validity. 

In FigJS]we compare different regimes of Eq.( !33l) with the exact numerically calculated IV- 
characteristic for forward bias. The low voltage part of all the IV-characteristics computed 
in Figs.([n]) and ([H]) for different values of the injection barrier can be understood in the 
following way. For voltages V < kT/e —Vbi the expression for the IV-characteristic, 
Eq. (l33l . can be approximated by 
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and shows a linear dependence of the current on the voltage. At zero voltage and current, 
electrons drift from the magnesium to the calcium electrode through the organic layer as 
a consequence of the difference in the work functions of the two electrodes. The blocking 
character of the right electrode restricts the effective electron density in the device as is seen 
from the exponential factor in Eq. (j34j) . 
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FIG. 8: IV- characteristics for an organic layer with thickness L = 100 nm contacted with a calcium 
electrode at x = —L/2 and a magnesium electrode at x = +L/2. The stars represent the exact 
numerically calculated IV- characteristic for forward bias. The solid line corresponds to the linear 
part approximated by Eq. (fff^U, the dashed line to the exponential current increase given by Eq. A35\) . 
the dotted line to the injection current described by Eq. and the dashed dotted line to a SCLC 
given by \j\ = (9/8) ee iJ,V 2 /L 3 . 



For voltages from the wide range kT/e <C V < —Vbi but not very close to — Vbi-, so 
that —e(V + Vbi) ^ kT, the IV-characteristic reveals a sharp exponential increase. Such a 
behavior can also be extracted by an approximation of equation fl33|) for voltages from this 
range (see Fig.([H])), 

The exponential increase stops when the voltage V approaches the value of the built-in 
voltage, —Vbi- At this point, the mean electric field in the device changes sign. 

Nevertheless, well above the built-in voltage, so that V + Vbi » kT/e, both the com- 
puted characteristic and Eq. (1331) show another exponential dependence of the current on the 
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voltage, due to the barrier lowering (see Fig.®): 

V + Vm. r ( A- lf F e(V + V B ,) \ 

In both equations, (J35l) and (|36l) . an exponential dependence on the voltage is observed. In 
the second case, however, this dependence is much weaker due to the factor l^ F /L e jf in the 
exponent. The latter exponential increase results from the nonvanishing width of the space 
charge zone in the injecting electrode. The appearance of this zone is a consequence of our 
self-consistent description of the device whole. 

After all, when the applied voltage overcomes the barrier at the injecting electrode, 
A - , the calcium contact can supply more charge carriers than the bulk of the organic 
semiconductor can transport. That is why all the current-voltage characteristics end up in 
the space charge limited regime with j ~ V 2 as is seen in Figs.([n]) and (|SJ). Particularly, 
when the calcium injection barrier is below the critical value of A cr u ~ 0.27 eV, the system 
shows a space charge limited current j ~ V 2 immediately after the built-in voltage is crossed, 
which can also be seen in Figj6j 



IV. CONCLUSIONS 



In this work, we have proposed a device model capable of describing a 
metal/insulator/metal device under injection limited as well as space charge limited condi- 
tions. The problem of defining the boundary conditions at the metal/insulator interfaces was 
solved by a self consistent treatment which fully includes the metal electrodes in a consistent 
one dimensional description of the device. In this treatment boundary conditions are defined 
far into the bulk of the metal electrodes where the respective media is well defined. The 
values for the electric field and the charge carrier densities at the interface can be calculated 
and depend on the condition of the considered system. 

We applied our model to organic semiconductors sandwiched between metal contacts, 
being interesting for optoelectronic applications. Though existing models include specific 
features of organic semiconductors, like energetic disorder or hopping transport, they de- 
scribe charge-carrier injection mostly in the two limiting cases of high or low injection barri- 
ers. In the former case charge- carrier injection models are based on the one-electron-picture 
within the framework of the Fowler- Nor dheim tunneling model or the Richardson-Schottky 
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model of thermionic injection , Q , where space charge effects are 

ignored completely. In the latter case the metal/insulator contact is assumed to be ohmic. 
Here, the contact can supply much more charge- carriers than can be transported through 
the bulk of the insulator and hence bulk properties alone define the characteristics of the 
system with no influence of the injection barrier JjJ, [2C|, |21[. As far as we know, the only 
model accounting consistently for both injecting barriers and space charge effects is that of 
Tut is et al. [la ]. The latter numerical treatment, however, does not allow analytical fitting 
of the current-voltage characteristics, making the physical reasons for the cross-over from 
barrier to space-charge dominated type of device behavior inaccessible. 

By contrast, the approach described here, offers the possibility to calculate IV- 
characteristics of metal/insulator/metal systems using experimentally accessible input pa- 
rameters. Our model is capable to predict IV-characteristics in forward and reverse bias 
and for all values of the injection barriers at each interface. The description includes the 
built-in voltage and field dependent effective injection barriers, being direct consequences 
of the self consistent approach. Calculated IV-characteristics can be divided in different 
regimes for which approximate solutions were derived. For low voltages a linear dependence 
of the current on the voltage is observed, followed by an exponential increase. This behavior 
is not due to intrinsic charge carriers, but due to diffusive charge carrier transport. Once 
the voltage exceeds the built-in voltage a weak exponential current increase is observed as a 
consequence of injection barrier lowering. Here, charge carrier injection dominates the IV- 
characteristic. At very high voltages, barrier lowering becomes so strong that more charge 
carriers are injected than can be transported through the bulk of the semiconductor. This 
leads to the SCL dependence, j oc V 2 , for all calculated IV-curves. The relative occurrence 
of these different regimes thereby depends strongly on the injection barriers at the contact 
between the insulator and the anode and the cathode metals, respectively. As a consequence, 
rectifying behavior is observed for strongly disparate contacts. 

A more realistic model incorporating charge-carrier trap states and Gaussian DOS distribu- 
tions of the organic semiconductor is currently developed. 
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